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Ever since the inception of gravitational-wave detectors, limits imposed by quantum mechanics to 
the detection of time-varying signals have been a subject of intense research and debate. Drawing 
insights from quantum information theory, quantum detection theory, and quantum measurement 
theory, here we prove lower error bounds for waveform detection via a quantum system, settling the 
long-standing problem. In the case of optomechanical force detection, we derive analytic expressions 
for the bounds in some cases of interest and discuss how the limits can be approached using quantum 
control techniques. 
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I. INTRODUCTION 

The study of quantum measurement has come a long 
way since the proposal of wavefunction collapse by 
Heisenberg and von Neumann, the philosophical debates 
by Bohr and Einstein, and the cat experiment hypothe- 
sized by Schrddinger. With more and more experimental 
demonstrations of bizarre quantum effects being realized 
in laboratories, many researchers have shifted their focus 
to the practical implications of quantum mechanics for 
precision measurements, such as gravitational-wave de- 
tection, optical interferometry, atomic clocks, and mag- 
netometry [l|-13]. Braginsky, Thorne, Caves, and others 
pioneered the application of quantum measurement the- 
ory to gravitational- wave detectors 0-Ql, while Holevo, 
Yuen, Helstrom, and others have developed a beautiful 
theory of quantum detection and estimation _8, 9] based 
on the more abstract notions of quantum states, effects, 
and operations [l3]- Although Holevo et aL's approach 
was able to produce rigorous proofs of quantum limits to 
various information processing tasks, so far it has been 
applied mainly to simple quantum systems with trivial 
dynamics measured destructively to extract static pa- 
rameters. Applying such an approach to gravitational- 
wave detection, or optomechanical force detection in gen- 
eral [llj, proved to be far trickier; the signal of interest 
there is time-varying (commonly called a waveform in 
engineering literature [f2]), the detector is a dynamical 
system, and the measurements are nondestructive and 
continuous [E-Q- Quantum limits to such detectors had 
been a subject of debate p^l - fTsj , with no definitive proof 
that any limit exists. In more recent years, the rapid 
progress in experimental quantum technology suggests 
that quantum effects are becoming relevant to metro- 
logical applications and has given the study of quantum 
limits a renewed impetus [11-0, [ll| . 

Generalizing the quantum Cramer-Rao bound first 



proposed by Helstrom [8|, Tsang, Wiseman, and Caves 
recently derived a quantum limit to waveform estima- 
tion jl6| , which represents the first step towards a rigor- 
ous treatment of quantum limits to a waveform sensor. 
That work assumes that one is interested in estimating 
an existing waveform accurately, so that the mean-square 
error is an appropriate error measure. The first goal 
of gravitational-wave detectors is not estimation, how- 
ever, but to detect the existence of gravitational waves, 
in which case the miss and false-alarm probabilities are 
the more relevant error measures [l2] and the existence 
of quantum limits remains an open problem. Here we 
settle this long-standing question by proving lower error 
bounds for the quantum waveform detection problem. To 
illustrate our results, we apply them to optomechanical 
force detection, demonstrating a fundamental trade-off 
between force detection performance and precision in de- 
tector position, and discuss how the limits can be ap- 
proached in some cases of interest using a quantum-noise 
cancellation (QNC) technique [T7l - [2^ and an appropriate 
optical receiver, such as the ones proposed by Kennedy 
and Dolinar [E, H^- Merging the continuous quantum 
measurement theory pioneered by Braginsky et al. and 
the quantum detection theory pioneered by Holevo et al. , 
these results are envisaged to play an influential role in 
quantum metrological techniques of the future. 



II. QUANTUM DETECTION OF A CLASSICAL 
WAVEFORM 



Let P[j/|'Ho] be the probability functional of an obser- 
vation process y{t) under the null hypothesis "Hqj and 



P[v\'Hi] = J DxP[x]P[y\x,Hi] 



(2.f) 
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be the probability functional under the alternative hy- 
pothesis Hi- x{t) is a classical waveform, P[x\ is its 
prior probability functional, and P[y\x,'Hi\ is the like- 
lihood functional under "Hi. To perform hypothesis test- 
ing given a record of y{t), one separates the observation 



space into two decision regions Tq and Ti, such that Hq 
is chosen if y falls in Tq and "Hi is chosen if y falls in Ti. 
The miss probability is defined as 



Poi = / DyP[y\Hi] 
and the false-alarm probability is 

Pio = / DyP[y\Ho]. 



(2.2) 



(2.3) 



Two popular decision strategies are the Bayes criterion, 
which minimizes the average error probability 



(2.4) 



given the prior hypothesis probabilities Pq and Pi — 1 — 
Pq , and the Neyman-Pearson criterion, which minimizes 
Poi for an allowable Pio, or vice versa [l2j . 

To introduce quantum mechanics to the problem, as- 
sume that x{t) perturbs the dynamics of a quantum sys- 
tem under "Hi and y(t) results from measurements of the 
system. Without any loss of generality, we model P[y|'Ho] 
and P[y\x,Hi] by considering a large enough Hilbert 
space, such that the initial quantum state \^) at time ti is 
pure, the evolution in the Schrodinger picture is unitary, 
and measurements arc modeled by a positive-operator- 
valued measure (POVM) E[y] at the the final time tf via 
the principle of deferred measurement [l^ [l^ HJ] : 

P[y\Ho]^tT{E[y]Uoitf,U)\ij){mUtf^U)}, (2.5) 
P[y\x,Hi]=tT[E[y]Uiitf,U)\i;){mUtf,U)}, (2.6) 

where only the unitaries Uq and Ui are assumed to differ 
and Ui depends on x. Assume further that 



Uo{tf,ti) = Texp 



Ui{tf,ti) = Texp 



dtHo{t) 



dtHi{x{t),t) 



Hi{x{t),t) ^ Hoit) + AH{xit),t), 



(2.7) 

(2.8) 
(2.9) 



where T denotes time-ordering and AH{x{t),t) is the 
Hamiltonian term responsible for the coupling of the 
waveform x{t) to the quantum detector. Figure [T] shows 
the quantum-circuit diagrams [1^ that depict the prob- 
lem. 

This setup can now be cast as a problem of quantum 
state discrimination between a pure state 



and a mixed state 



Pi 



DxP[x]Ui\'ip}{^/j\Ul. (2.11) 
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FIG. 1. Quantum-circuit diagrams for the waveform detection 
problem. The quantum system is modeled as a pure state \tp} 
with unitary evolution {Uo or Ui) under each hypothesis ("Ho 
or Hi) in a large enough Hilbert space for a given classical 
waveform x{t), which perturbs the evolution under Hi. If 
x{t) is stochastic, the final quantum state under Hi is mixed. 
Measurements are modeled as a positive-operator- valued mea- 
sure (POVM) E[y] at the final time through the principle of 
deferred measurement. 



Let l^j) be a purification of pj in a larger Hilbert space 
(gjHs, such that pj ^ trs and tT{E[y]pj} = 

tr[(£:[y] ® 1b)|*j)(*j|], where 1b denotes the identity 
operator with respect to Hb. The average error proba- 
bility is thus lower-bounded by Q: 



I 



Pe > 2 (l - v/l-4PoPi|('J'o|*i)pj , 
which is valid for any purification. Hence 



Pe> X I- /l-4PoPi max |(vI/o|*i)r 
2 \ V li-o),!*!) 



l-4PoPiP , 



(2.12) 

(2.13) 
(2.14) 



where F is the quantum fidelity by Uhlmann's theorem 



P(po, Pi) = (^tr ^/plpo^/pl 
As Pq is pure, the fidelity is given by 

P=(V^|[/oViC/o|^)=E(P,), 



(2.15) 

(2.16) 
(2.17) 
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where we have defined classical and quantum averages by 

(2.18) 
(2.19) 



li-) = J DxP[x]i-), 



By similar arguments, a quantum bound on the miss 
probability Pqi for a given allowable false-alarm probabil- 
ity Pin can be derived from the bound for the pure-state 
case 



0, 



Pio > F. 
(2.20) 



Note that the latter bound is equally valid if we inter- 
change Pqi and Pio; for example, fixing Pqi = means 
Pio > P. Equations (|2.14p and (|2.20p are valid for any 
POVM and achievable if x(t) is known a priori, such that 
both po and pi are pure Q. 

In terms of related prior work at this point, Ou [26j and 
Paris [131 studied quantum limits to interferometry in the 
context of detection, while Childs et al. [1^ , Acm et al. 
[l^HO], and D'Ariano et al. [3l| also studied unitary or 
channel discrimination, but all of them did not consider 
time-dependent Hamiltonians, which are the subject of 
interest here. 

A key step towards simplifying Eq. (|2.17p is to recog- 
nize that 



ul{tf,U)Ui{tf,U)=T exp 



dtAHo{x{t),t) 



(2.21) 



where 



AHoix{t),t) = W^{t,U)AH{x{t),t)Uo{t,U) (2.22) 

is AH in the interaction picture Q . In general, Eq. (I2.17P 
can then be expanded in a Dyson series and evaluated 
using perturbation theory [32|. To derive analytic ex- 
pressions, however, we shall be more specific about the 
Hamiltonians and the initial quantum state. 



III. FORCE DETECTION WITH A LINEAR 
GAUSSIAN SYSTEM 

Assume that a; is a force on a quantum object with 
position operator q, so that 

AH = -qx, (3.1) 

and the conditional fidelity P^ becomes 



Texp 



dtqo{t)x{t) 



(3.2) 



with qo(t) obeying equations of motion under the null hy- 
pothesis Hq in the interaction picture. The (•) expression 



in Eq. (|3.2I) is a noncommutative version of the character- 
istic functional f33[. To simplify it, assume further that 
Hq consists of terms at most quadratic with respect to 
canonical position or momentum operators, such that the 
equations of motion are linear and qo (t) depends linearly 
on the initial-time canonical operators. Let Z{t) be a col- 
umn vector of canonical position/momentum operators, 
including qo{t), that obey the equation of motion 



dZ{t) 
dt 



= G{t)Z{t) + J{t) 



(3.3) 



under hypothesis "Hq, where G{t) is a drift matrix and 
J(t) is a source vector, both consisting of real numbers. 
qo(t) can then be written as 

qo{t)^Vq{t,ti)Z{t,)+ dTVgit,T)J{T), (3.4) 

where Vq(t,ti) is a row vector and a function of G(t). 
This gives 



1 



dtqo{t)x{t) = KJZ{ti) -t- 0, 



dtx{t)Vqit,t,), 



dtx{t) / dTVq{t,T)J{T). 
ti Jti 



With Fx now given by 

Fx = I (Texp [zK^Z(t,)+i0])f , 
the time-ordering operator becomes redundant: 

Fx = |(V'|exp[zK^Z(<,)] IV^)!'. 



(3.5) 
(3.6) 
(3.7) 

(3.8) 
(3.9) 



This expression can be simplified using the Wigner rep- 
resentation W{z, ti) of It/)), which has the following prop- 
erty [ail: 



I exp [iK^Z{ti)\ 



= I dzW{z,ti) exp{iK z), 

(3.10) 



where z is a column vector of phase-space variables. As- 
suming further that W{z,ti) is Gaussian with mean vec- 
tor z and covariance matrix S, we obtain an analytic 
expression for Fx : 



Fx 



dzW{z, ti) exp(zK^z) 
exp (— k^Ek) 



(3.11) 
(3.12) 



exp 



dt 



dt'x{t)Y.q{t,t')x{t') 



^q{t,t') = Vq{t,ti)^V;j {t' ,ti). 



(3.13) 
(3.14) 
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The covariance matrix is given by the Weyl-ordered sec- 
ond moment; 



^jk = 2 {Zj{ti)Zk{ti) + Zk{ti)Zj{ti)) 
-{Z,{U)) {ZkiU)). 



(3.15) 



Hence 



^q{t,t') = - {qo{t)qo{t') + qo{t')qo{t)) - {qo{t)) {qo{t')) . 

(3.16) 

It is interesting to note that the expression given by 
— In in Eq . p.l3p coincides with the one proposed 
in Refs. |35| as an upper quantum hmit on the force- 
sensing signal-to-noise ratio, and AYiq{t,t')/'h? is equal to 
the quantum Fisher information in the quantum Cramer- 
Rao bound for waveform estimation [16]. The relation 
of this expression to the fidelity and the detection error 
bounds is a novel result here, however. 

If the statistics of qo{t) can be approximated as sta- 
tionary; viz.. 



/oo 
-^Sq{uj)eM-Mt-t')] 
-oo 2vr 



Fx becomes 

F^ = exp 

x{uj) = / dtx{t) exp(jcji) 



For example, if 

is a sinusoid. 
Fx « exp 



x{t) = Xcoii{VLt + I 



T = tf-U. 



(3.17) 

(3.18) 
(3.19) 

(3.20) 

(3.21) 



These expressions for the fidelity suggest that, for a given 
x{t), there is a fundamental trade-off between force de- 
tection performance and precision in detector position. 



IV. OPTOMECHANICS 

Suppose now that the mechanical object is a moving 
mirror of an optical cavity probed by a continuous-wave 
optical beam, the phase of which is modulated by the ob- 
ject position and the intensity of which exerts a measure- 
ment backaction via radiation pressure on the object, as 
depicted in Fig. [21 This setup provides a basic and often 
sufficient model for more complex optomechanical force 
detectors. Let the output field operator under hypothesis 
be 




FIG. 2. A cavity optomechanical force detector. An optical 
cavity with a moving mirror is pumped on-resonance with an 
input field Ain, and the output field ylout is measured to infer 
whether a force x{t) has perturbed the motion of the mirror. 



where 

a{t) * b{t) 
denotes convolution. 



dw 
2^ 



dTa{t - T)b{T) 



Kn{uj) ejq){—iu!t) 



is an impulse-response function with 

iuj + J 



K2H = 



—iuj -f 7 
2a;o 1 

L —iuj + 7 



(4.2) 
(4.3) 

(4.4) 
(4.5) 



in the frequency domain, A is the input mean field, wq 
is the optical carrier frequency, L is the cavity length, 
and 7 is the optical cavity decay rate [l3]. qj{t) is the 
position operator under each hypothesis, which can be 
written as [TtI 

qo{t)^K:i{t)*hK2{t)*m, (4.6) 
gi(i) « if3(t) * [nK2{t) * m + <t)l (4.7) 

where K'i{t) is another impulse response function that 
transfers a force to the position, 



(4.8) 



is the backaction noise, and the transient solutions are 
assumed to have decayed to zero. Defining 



Ki{t) = K:,{t)*K2{t), 
such that the position power spectral density is 

we obtain 

do; 



A 



ontj 



{t) « Ki{t) * Ain(i) + tAK2it) * q,it), (4.1) 



Fx 



exp 



2tt 



S^{u)\K4u)x{lu)\' 



(4.9) 
(4.10) 

(4.11) 
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The backaction noise ^ that appears in the output 
field, in addition to the shot noise in Ai^, can limit the 
detection performance at the so-called standard quan- 
tum limit 0-0, [iBl- This does not seem to agree with 
the fundamental quantum limits in terms of Eq. (j4.1ip . 
which suggest that increased fluctuations in qo{t) due to 
^(t) can improve the detection. Fortunately, it is now 
known that the backaction noise can be removed from 
the output field [H M, M^^- One method, called 
quantum- noise cancellation (QNC), involves passing the 
optical beam through another quantum system that has 
the effective dynamics of an optomechanical system with 
negative mass [13, [H, [20l - [2^ . With the backaction noise 
removed, the output fields become 

^outO 

(<) «i^i(i)*A„(<), (4.12) 
Aontiit) « Ki{t) * Ai„(<) + iAK2it) * Kait) * 

(4.13) 

If the phase quadrature of Ao^tjit) is measured by ho- 
modyne detection, the outputs can be written as 



yoit) « r]{t), 

yi{t) ^ r]{t) + K2{t) * Ksit) * x{t) 



2i\A\ 



(4.14) 
(4.15) 

2 ,A*Ki{t) * AAi„(t) - AKlit) * AAj-Jt)]. 

(4.16) 



The power spectral densities of ^(i) and 77 (i) satisfy an 
uncertainty relation [5] : 



(4.17) 



The detection problem described by Eqs. (|4.14l) and 
(|4.15p becomes a classical one with additive Gaussian 
noise, a scenario that has been studied extensively in 
gravitational- wave detection [s^, HI] . 



V. ERROR BOUNDS FOR DETERMINISTIC 
WAVEFORM DETECTION 

Suppose that x{t) is known a priori. It is then 
well known that the error probabilities for the detec- 
tion problem described by Eqs. (|4.14p and (14. 15^ using a 
likelihood-ratio test are [12] 



where 



-PlO,hom = 2 ^^''^^ V ^ 4^ 

-Poi.hom = ^ erfc fcr - ^ 



2 f°° 

erfcu= / dvexp{—v ), 

Ju 



(5.1) 
(5.2) 

(5.3) 



A is the threshold in the likelihood-ratio test, which can 
be adjusted according to the desired criterion, and ct is a 
signal-to-noise ratio given by 



1 



dui \K4{uj)x{lu)\^ 



2tt 



S',,(w) 



(5.4) 



for a long observation time relative to the duration of 
x{t) plus the decay time of K/^{t). To compare homo- 
dyne detection with the quantum limits, suppose that 
the duration of x{t) is long and increases at least lin- 
early with T , so that we can define an error exponent as 
the asymptotic decay rate of an error probability in the 
long-time limit. For simplicity, we consider here only the 
exponent of the higher error probability: 



lim — lnmax{Plo,Pol}• 
T— >oo 1 



(5.5) 



Although this asymptotic limit may not be relevant to 
gravitational- wave detectors in the near future, the error 
probabilities for which are anticipated to remain high, 
we focus on this limit to obtain simple analytic results, 
which allow us to gain useful insight into the fundamental 
physics. More precise calculations of error probabilities 
are more tedious but should be straightforward following 
the theory outlined here. 

For homodyne detection, the error exponent is 



horn 



1 

8T 



duj \K4^(lu)x{lu)Y 



27r 



Sriiuj) 



(5.6) 



The quantum limit, on the other hand, is 

- lim ;^lnmax{Pio,Poi} 
< — lim — In max P„ < — lim — \nF 

- T~ioo T Pa T^oo T 



which gives 



Tf = - 
T 



'^S^{u:)\K,{u:)x{u:)\^. 
Ztt 



(5.7) 



(5.8) 



Using the uncertainty relation between and in 
Eq. (|4.17p . it can be seen that 



(5.9) 



that is, the homodyne error exponent is at most half 
the optimal value. This fact is well known in the con- 
text of coherent-state discrimination j8|, i23, [s^ |4^ . The 
suboptimality of homodyne detection here should be con- 
trasted with the conclusion of Ref. , which states that 
homodyne detection together with QNC are sufficient to 
achieve the quantum limit for the task of waveform esti- 
mation. 

To see how one can get closer to the quantum limits, 
let's go back to Eqs. (|4.12l) and (|4.13p . Observe that, if 
the input field is in a coherent state, the output field is 



also in a coherent state (in the Schrodinger picture) un- 
der each hypothesis. This means that existing results for 
coherent-state discrimination can be used to construct 
an optimal receiver. The Kennedy receiver, for exam- 
ple, displaces the output field so that it becomes vacuum 
under 'Hn and then detects the presence of any output 
photon Q. Any detected photon means that Hi must 
be true. Deciding on Hq if no photon is detected and Hi 
otherwise, the false-alarm probability Piq is zero, while 
the miss probability is the probability of detecting no 
photon given "Hi, or 




Photon detector 



out 



displace by —A. 



01, Ken = exp 



dt\AK2{t)*K3{t)*x{t)\' 



(5.10) 



For a long observation time with = \A\'^ for a coherent 

state, 



-Poi^Kcn ~ exp 



OO J 



F, (5.11) 



which makes the Kennedy receiver optimal under the 
Neyman-Pearson criterion in the case of Pio = ac- 
cording to Eq. (|2.20l) and also achieve the optimal error 
exponent: 

Tkoii = - lim T^lnPoi.Kcn = Tf. (5.12) 

The Kennedy receiver can be integrated with the QNC 
setup; an example is shown in Fig. [3] The Dolinar re- 
ceiver, which updates the displacement field continuously 
according to the measurement record, can further im- 
prove the average error probability slightly to saturate 
the lower limit given by Eq. (|2.14l) 0, [2I] . Other more 
recently proposed receivers may also be used here to beat 
the homodyne limit [s^, H^] . 



VI. ERROR BOUNDS FOR STOCHASTIC 
WAVEFORM DETECTION 

Consider now a stochastic x{t), which should be rele- 
vant to the detection of stochastic backgrounds of gravi- 
tational waves l4l|. Since is Gaussian, 



F = / DxP[x\ exp 



-■^ / dtdt'x{t)Y.q{t,t')x{t') 



(6.1) 



can be computed analytically if the prior P[x\ is also 
Gaussian. Here we shall use a discrete-time approach and 
take the continuous limit at the end of our calculations. 
If x{t) is a zero- mean Gaussian process with covariance 

^^{t,t') = nx{tW)l (6.2) 



FIG. 3. An integrated QNC-Kennedy receiver. The output 
field Aout from the optomechanical force detector in Fig. [2] is 
displaced by —A and then passed through an optical setup 
that removes the measurement backaction noise. The dash 
arrows represent a red-detuned optical cavity mode that mim- 
ics a negative-mass oscillator and interacts with the optical 
probe field via a beam splitter (BS) and a two-mode optical 
parametric amplifier (OPA). Details of how this setup works 
can be found in Refs. [13, [lal- If the field Ain is in a coherent 
state, the final output field should be in a vacuum state under 
the null hypothesis Ho- Any photon detected at the output 
indicates that T-ii must be true. 



it can be discretized as 

x= {xn,. . . ,xn-iY , 
DxP[x\ w dxQ . . . dxN^i 



V(2^)^dctS, 



1 



X exp ( --x'^T,^^x 



= E{xx'^). 



(6.3) 

(6.4) 
(6.5) 



The fidelity then becomes a finite-dimensional Gaussian 
integral: 



F fa dxo . . . dxN-i 



1 



v/(27r)^ det 

/ 1 T -1 T 
\ 2 tf ^ 

/det(S^i + 2(5i2s,/n2)-i 



det 

det ( / + ^^EqEa; 



exp 



-1/2 



- - tr In I / -I- -^-HqY^x 



1, . f. 25^ 
2 



where A^^ are the eigenvalues of the matrix 



(6.6) 

(6.7) 

(6.8) 
(6.9) 
(6.10) 

(6.11) 
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If T,q{t,t') and ^x{t,t') are both stationary; viz., 

^git,t') = a,it-t'), (6.12) 
^At,t')^ax{t~t'), (6.13) 

they can be modeled as circulant matrices in discrete 
time, so that C is also circulant, with eigenvalues given 
by the discrete Fourier transform of a row or column 
vector of the matrix. Taking the continuous-time limit 
using 



at the output. Although such assumptions are highly 
idealistic, especially for current gravitational-wave detec- 
tors, the ideal scenario shows that the quantum bounds 
proposed here are in principle achievable using known 
optics technology. Optimal discrimination of squeezed or 
other Gaussian states remains a topic of current research 
[3 - |46j and may be useful for future gravitational-wave 
detectors that use squeezed light j47,] . Generalization of 
the results here to multi- waveform discrimination should 
also be useful for gravitational-wave astronomy [s^, HI] 
and may be done by following Refs. ^, A9] . 



(6.14) 



we get 



F 



1 duj^ 
— / — m 

2 ./ _ 27r 



1 + ■^Sq{uj)Sa;{uj) 



dt<Tx{t) exp(iaji). 



(6.15) 
(6.16) 

(6.17) 

(6.18) 



This fidelity expression can then be used in the detection 
error bounds. 

For homodyne detection, the error exponent is more 
complicated for stochastic waveform detection and given 
by the Chernoff distance [42|, |43|]: 



J- horn 

1 r diJ l + {l-s)\Kiiuj)\^Sxiuj)/S^{uj) 



(6.19) 



The performance of homodyne detection relative to the 
quantum limits then depends on the specific form of 
\K4{uj)\'^Sx{oj)- The Kennedy receiver, on the other 
hand, is still applicable here, as the output is still a co- 
herent state under Hq . The false-alarm probability is still 
zero, and the miss probability is now 



-Poi,Kcn ~ Eexp 



(6.20) 



which means that the Kennedy receiver remains opti- 
mal, both in terms of the Neyman-Pearson criterion in 
the case of Pio = and the error exponent. Whether 
other receivers can do even better and saturate the other 
quantum bounds is a more difficult question, as the out- 
put field under "Hi is now in a mixed state and the fidelity 
lower bounds may not be achievable. 

The use of Kennedy or Dolinar receivers assumes co- 
herent states at the output, which is the case only if the 
backaction noise cancellation is complete and quantum 
shot noise in the input beam is the only source of noise 



VII. OUTLOOK 

Now that quantum limits to waveform detection have 
been discovered, the natural next question to ask is how 
they can be approached in practice. In the case of op- 
tomechanical force detection, the requirements are quan- 
tum shot noise as the only source of noise at the output 
and an appropriate receiver, such as the Kennedy re- 
ceiver. A proof-of-concept experimental demonstration 
of waveform detection approaching the shot-noise limits 
based on Eq. (|4.11l) should be well within reach of current 
quantum optics technology. To demonstrate the trade- 
off between force detection performance and detector lo- 
calization suggested by Eq. p.l8p with optomechanics, 
however, can be much more challenging, as it would re- 
quire quantum backaction noise to dominate the detector 
position fluctuation but become negligible in the output 
via QNC. A more promising candidate for this demon- 
stration is atomic spin ensembles, with which backaction- 
noise-canceled magnetometry has already been realized 
[23 |. The likelihood-ratio formulas derived in Ref. [s^l 
should be used in practice instead of the ideal-case deci- 
sion rules discussed here to account for any excess noise. 

In terms of potential further theoretical work, it should 
be useful to generalize beyond the assumptions of scalar 
waveform, linear Gaussian systems, optical coherent 
states, stationary processes, and long observation time 
used here. The fidelity expressions derived here r nay also 
be useful for the study of waveform estimation [la], ei- 
ther as an alternative way of deriving the quan tum Fisher 
information via a Taylor-series expansion_|5l| or used di- 
rectly in the quantum Ziv-Zakai bound (52| . 

From a more conceptual point of view, this study, to- 
gether with the earlier work on waveform estimation ilQ] , 
shows that the concepts of states, effects, and opera- 
tions fade into background when dealing with dynami- 
cal quantum information systems, and multi-time quan- 
tum statistics, through the use of Heisenberg or interac- 
tion picture, take the center stage. It may be interesting 
to explore whether this perspective has any relevance to 
other dynamical quantum information systems, such as 
quantum computers [24*1, and the study of quantum cor- 
relations |53i] . 
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